EXISTENCE AND BLOW-UP BEHAVIOR FOR SOLUTIONS 
OF THE GENERALIZED JANG EQUATION 



QING HAN AND MARCUS KHURI 

Abstract. The generalized Jang equation was introduced in an attempt to prove the Penrose in- 
equality in the setting of general initial data for the Einstein equations. In this paper we give an 
extensive study of this equation, proving existence, regularity, and blow-up results. In particular, 
precise asymptotics for the blow-up behavior are given, and it is shown that blow-up solutions are 
not unique. 



1. Introduction 

In 1978 the physicist P.-S. Jang [12J introduced a quasilinear elliptic equation in connection with 
the positive energy conjecture of General Relativity. Since then it has had numerous applications 
from quasilocal mass to the existence of black holes. An excellent survey of these applications may 
be found in [2j. However, Jang's equation has been shown to be unsuitable for an application to the 
Penrose inequality |14j . and for this reason a generalization of this equation has been proposed in [1] 
and [5]. A further application may be found in [13j . The purpose of this paper is to give a complete 
analysis of the generalized Jang equation, in a setting suitable for the Penrose inequality. 

Consider an initial data set (M, g, k) for the Einstein equations. This consists of a Riemannian 
3-manifold M with metric g, and a symmetric 2-tensor k, which satisfy the constraint equations: 

2^ = R+{Trkf-\k\\ 
J = div(k + (Trk)g). 

Here \x and J are the energy and momentum densities of the matter fields, respectively, and R is the 
scalar curvature of g. If all measured energy densities are nonnegative then [i > \J\, which will be 
referred to as the dominant energy condition. This condition may be viewed as a quasi nonnegativity 
of the scalar curvature. In fact in the time symmetric (k = 0) or maximal case (Trk = 0), one does 
have R > 0, and this condition plays a central role in the proof of several results such as the positive 
energy theorem [16], |18| and the Penrose inequality [3], [11) . Moreover, the primary difficulty in 
establishing these results for general initial data is the lack of nonnegative scalar curvature. Thus, 
in the general case, one is motivated to deform the initial data in an appropriate way (depending 
on the problem at hand) such that the dominant energy condition yields nonnegativity of the scalar 
curvature for the deformed metric. The deformation chosen by Jang [12] is given by ~g = g + df 2 , for 
some function / defined on M. Notice that g is the induced metric on the graph E = {t = f(x)} in 
the product 4-manifold (M xR,g + dt 2 ). When viewed in this way, the optimal Jang deformation is 
given by a hypersurface E C (M xR,g + dt 2 ) which has scalar curvature that is as nonnegative as 
possible. It turns out that this optimal deformation arises as the hypersurface which satisfies Jang's 



The first author acknowledges the support of NSF Grant DMS-1105321. The second author acknowledges the 
support of NSF Grant DMS-1007156 and a Sloan Research Fellowship. 

1 



2 



HAN AND KHURI 



equation 

(1.2) # E - Tr^k = 0, 

where H% is mean curvature and Tr-^k denotes the trace of k (extended trivially to the 4-manifold) 
over S. While this deformation does not necessarily yield nonnegative scalar curvature, it is sufficient 
to make a further conformal deformation to zero scalar curvature, as was carried out by Schoen and 
Yau [17] in their proof of the positive energy theorem. Jang's approach also beautifully handles the 
case of equality for the positive energy theorem. In the case of equality, the Jang surface (£,g) is 
isometric to Euclidean space (M 3 ,<5), and hence g = 5 — df 2 . It follows that the map x h-» (x,f(x)) 
yields an isometric embedding of (M,g) into the Minkowski spacetime M 4 ; it can also be shown that 
the second fundamental form of this embedding agrees with k, as desired. In some sense, Jang's 
procedure seems to be tailor made for the positive energy theorem. In fact it is so well calibrated 
to the positive energy theorem, that this method is rendered inapplicable to the Penrose inequality. 
The reason for this is as follows. As explained in [?] and [5], in the case of equality for the Penrose 
inequality, the Jang metric g should coincide with g sc , the induced metric on the t = slice of the 
Schwarzschild spacetime §C 4 . It follows that g = g sc — df 2 , and the map 14 (x, f(x)) cannot yield 
an embedding of (M,g) into SC 4 , since the Schwarzschild metric has the warped product structure 
—cj} 2 c dt 2 + g sc . However, this observation leads to a natural modification of the Jang approach from 
which the generalized Jang equation arises. 

We now search for a hypersurface E, given by a graph t = f(x), inside the warped product space 
(M xR,g + 4> 2 dt 2 ), where <\> is a nonnegative function defined on M. For certain applications the 
choice of <p may depend on /, however in this paper we will assume that cfi is fixed independent of /. 
The goal or motive which leads to the generalized Jang equation is the same as in the classical case. 
That is, we search for a hypersurface which has the most positive scalar curvature that is possible. 
In order to have any chance of obtaining a positivity property for the scalar curvature, we would like 
the Jang surface E to satisfy an equation with the same structure as in (|1.2|) . namely 

(1.3) #s - TrxK = 0, 

where again Hj] is mean curvature and Tr^K denotes the trace of K over E. Here, however, K 
represents a nontrivial extension of the initial data k to all of M x M. In (|1.2|) k is extended trivially 
in that k{dt,-) = 0, but it turns out that the trivial extension is not appropriate for a Jang surface 
inside the warped product metric. This is due to the fact that for applications it is desirable for 
solutions of the generalized Jang equation to blow-up at apparent horizons, and as is shown in [3], 
this is possible when K is extended as follows: 



(1.4) 



K{d xl , d xj ) = K(d x3 ,d x i) = k(d x i ,d xJ ) for 1 < *,j < 3, 
K(d xi , d t ) = K(d t , d x i) = for 1 < i < 3, 



K(d t ,dt 



2 g(V/,V0) 



where x l , i = 1,2,3, are local coordinates on M. Moreover this particular extension also yields an 
optimal positivity property for the scalar curvature of solutions to ()1.3|) . Namely, a long calculation 
(|I]> [5]) gives the following formula for the scalar curvature of S satisfying equation (|1.3p with 
extension (|1.4[) : 



(1.5) R = 2(fi- J(w)) + \A- K\y,\ 2 + 2\q\ z - 20" 1 div(<f>q). 
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Here g = g + 4> 2 df 2 and A are the induced metric and second fundamental form of X, respectively, 
is the restriction to X of the extended tensor K, div is the divergence operator with respect to 
~g, and q and w are 1-forms given by 



with f 3 = g VJ fi. If the dominant energy condition is satisfied, then all terms appearing on the 
right-hand side of (|1.5p are nonnegative, except possibly the last term. However the last term has 
a special structure, and in many applications it is clear that a specific choice of </> will allow one to 
'integrate away' this divergence term, so that in effect the scalar curvature is weakly nonnegative 
(that is, nonnegative when integrated against certain functions). 

When the tensor k is extended according to (jl.4p , we will refer to equation (j 1 . 3|) as the generalized 
Jang equation, and the solution £ = {t = f(x)} will be called the Jang surface. In local coordinates 
the generalized Jang equation takes the following form: 



'-.CP \ / oV-,/ • o,.!) ■ O jJ) 

i + ^lv/lV^ y/i + <p\vff 

This is a quasilinear elliptic equation, which degenerates when / blows-up or <f> = 0. When <p = 1 
this reduces to the classical Jang equation studied by Schoen and Yau [17] , 

In this paper we will always assume that the initial data are asymptotically flat (with one end), 
so that at spatial infinity the metric and extrinsic curvature satisfy the following fall-off conditions 

\d m { gij -5 ij )\ = 0(\x\~ m - 1 ), \d m k ij \=0(\x\- m - 2 ), to = 0,1, 2, as \x\ oo. 

Moreover if the initial data has boundary, it will be assumed to consist of an outermost apparent 
horizon. To explain this more precisely, recall that the strength of the gravitational field in the 
vicinity of a 2-surface S C M may be measured by the null expansions 

(1.7) 9± := H s ± Tr s k, 

where Hs is the mean curvature with respect to the unit outward normal (pointing towards spatial 
infinity). The null expansions measure the rate of change of area for a shell of light emitted by the 
surface in the outward future direction (#+), and outward past direction (#-). Thus the gravitational 
field is interpreted as being strong near S if 6 + < or < 0, in which case S is referred to as a future 
(past) trapped surface. Future (past) apparent horizons arise as boundaries of future (past) trapped 
regions and satisfy the equation 6 + = (9- = 0). In the setting of the initial data set formulation 
of the Penrose inequality, apparent horizons take the place of event horizons, in that the area of the 
event horizon is replaced by the least area required to enclose the outermost apparent horizon. Here, 
an outermost future (past) apparent horizon refers to a future (past) apparent horizon outside of 
which there is no other apparent horizon; such a horizon may have several components, each having 
spherical topology ([7], [8]). In this paper we will refer to the union of the outermost future apparent 
horizon and outermost past apparent horizon as the outermost apparent horizon, and we will assume 
that past and future horizon components do not intersect. For the sake of applications it is desirable 
for solutions of the generalized Jang equation to blow-up to +oo (— oo) in the form of a cylinder over 
the future (past) components of an outermost apparent horizon, and to vanish at spatial infinity. 
The fall-off rate for the solution of the generalized Jang equation depends on the asymptotics of the 
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warping factor, however here we will always assume that 

(1.8) <j>(x) = 1 + + ( ^ ) as Id -)• oo 
for some constant C, which yields 

(1.9) |V m /|(x) = 0{\x\^- m ) as \x\-too, m = 0,1,2. 



These asymptotics ensure that the ADM energy of the Jang surface agrees with that of the initial 
data. 

Our first theorem generalizes the basic existence result of Schoen and Yau |17j . 

Theorem 1.1. Suppose that (M,g,k) is a smooth, complete, asymptotically flat initial data set, and 
that (ft is smooth, strictly positive, and satisfies (|1.8p . Then there exist disjoint open sets ri + , f2_ C M , 
and a smooth function f : M — (0+ U f2_) — > R satisfying the generalized Jang equation (|1.6|) as well 
as (|1.9j) . Furthermore dfl+ (8Q-) is a future (past) apparent horizon with f(x) — > ±oo as x —> d$l±. 
More precisely graph(/) is asymptotic to the cylinders <9f2 + x M + and x R_ . 

This result shows that solutions of the generalized Jang equation behave in a similar fashion to 
those of the classical Jang equation, as long as the warping factor eft is positive. If, however, cf) is 
allowed to vanish, which is required in some applications, then distinctly different behavior can arise. 

In Theorem II. 1\ the solution of the generalized Jang equation blows- up at apparent horizons, 
although which apparent horizons the solution blows up at is not indicated. For applications, how- 
ever, it is desirable to prescribe blow-up at the outermost apparent horizon. In fact, the existence of 
blow-up solutions for the classical Jang equation has been proven by Metzger in [13] (see also [6]), 
where he also establishes an upper bound for these solutions. In this paper we prove the existence 
of blow-up solutions for the generalized Jang equation, and establish upper and lower bounds for 
the solutions. In order to estimate the rate of blow-up for the solution at the horizon we will need 
to know the asymptotics of the null expansions as well. Set t(x) = dist(x,dM), and let S T denote 
level sets of the geodesic flow emanating from dM, that is, each point of S T is of distance r from 
the boundary. Furthermore decompose dM = d+M U d-M, where 9+M (d-M) denotes the future 
(past) apparent horizon components. We then stipulate that near d±M 

(1.10) c-V < 9±(S T ) < ct 1 , 

for some constants I, c > 0. 

Theorem 1.2. Suppose that (M,g,k) is a smooth, asymptotically flat initial data set, with outermost 
apparent horizon boundary dM = d+M U cLM . Suppose further that (ft is smooth, strictly positive, 
and satisfies (jl.8p . Ifl>l and (jl.lOp is valid, then there exists a smooth solution f of the generalized 
Jang equation (jl.6p . satisfying (jl.9p . and such that f{x) — > ±oo as x — >■ d±M . More precisely, 
graph(/) is asymptotic to the cylinders d+M x R + and d-M x R_, and in a neighborhood of d±M 



a -1 r"V + P~ l < ±f < ar' 1 -^ + (3 if I > 1, 
(1.11) - J ~ 

-orMogT + /T 1 < ±f < -alogr + /3 if 1 = 1, 

for some positive constants a and f3. 



Although this theorem gives the desired blow-up at the outermost apparent horizon, it is not 
sufficient for application to the Penrose inequality since the warping factor is not allowed to vanish. 



THE GENERALIZED JANG EQUATION 



5 



In the next result we remove this restriction, assuming that the warping factor vanishes in a controlled 
way at the horizon. Namely, in a neighborhood of dM we impose the following structure condition, 

(1.12) 4>{x) = T b (x)0(x) 

for some smooth (up to the boundary), strictly positive function <ft, with b > 0. The bounds in 
(jl.lOP are also quite restrictive, however these will also be relaxed in the next theorem for one of the 
asymptotics. In particular we will only require that near d±M, 

(1.13) |0±(S T )| < ct 1 

for some constants l,c > 0. Notice that when the initial data are smooth up to the boundary, as 
will always be the case in this paper, I > 1. Whether or not / actually blows- up and approximates a 
cylinder over the horizon is highly dependent on the relationship between the vanishing rates of 6± 
and 4>, that is, the relation between I and b. 

Theorem 1.3. Suppose that (M,g,k) is a smooth, asymptotically flat initial data set, with outermost 
apparent horizon boundary dM = d+M U d-M . Suppose further that (p is smooth, strictly positive 
away from dM , and satisfies (|1.8p and (I1.12p for some b > 0. 

(1) If - l -=± <b< l -±± and (fTTTOl) is valid for some I > 1, then there exists a smooth solution f of 
the generalized Jang equation (jl.6p . satisfying (jl.9p . and such that f(x) — > ±oo as x — > d±M . More 
precisely, in a neighborhood of d±M 

a- 1 T- b -^+r 1 <±f<aT- b -^+p if - l -^<b< l -±±, 

(1.14) 2 2 
-a^logr + zT 1 < ±f < -alogr + /3 if b = — , 

for some positive constants a and (5. 

(2) If | < b < and (fTT3|) is valid for some I > 1, then there exists a smooth solution f of 
the generalized Jang equation (jl.6p . satisfying (|1.9p . and such that f(x) — > ±oo as x — > d±M . More 
precisely, in a neighborhood of d±M 



(1.15) 



-a~ l \ogT + ft' 1 <±f <-a\ogT + ft if b = -, 



for some positive constants a and f3. 

Remark 1.4. Note that the hypothesis b > — is trivially satisfied under the primary assumptions 
of this paper, that is b > and I > 1. However it is believed that this theorem, and others having 
this hypothesis, continue to hold even when these primary assumptions are relaxed. For this reason, 
the inequality b > — -^ ls included in the statement of such results. 

We note that there are at least two solutions with different asymptotics when ^ < b < Mp. The 
estimates (I1.14p and (]1.15p are established by constructing appropriate sub and super solutions, 
however these estimates may not hold when the particular inequalities between I and b are not 
satisfied. Furthermore, we exhibit by example, a solution which does not have cylindrical asymptotics 
near the outermost apparent horizon when b = (I + l)/2. 

Theorem 11.31 may be considered the primary result of this paper, as it establishes the existence 
of solutions to the generalized Jang equation which possess appropriate behavior for application to 
the Penrose inequality. In particular, for an appropriate choice of 0, the blow-up rates (|1.15p show 
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that the Jang surface is a manifold with boundary, and that the boundary is a minimal surface. 
This allows techniques from the time symmetric proof of the Penrose inequality to be applied, and 
is different from the behavior of the blow-up solutions of the classical Jang equation which possess 
an infinitely long neck at the horizon. 

This paper is organized as follows. In Section [2] we derive local a priori estimates, and in Section [3] 
an important Harnack inequality is proved. These two results are then combined with further global 
estimates to establish Theorem 11.11 in Section HI In Section [5] we produce the blow-up solutions of 
Theorem 11.21 and construct appropriate super solutions. Lastly, the more difficult sub solutions are 
constructed in Section [6] and Theorem 11.31 is proven. 

2. A Priori Estimates for the Generalized Jang Equation 

In this section we begin the proof of existence for the generalized Jang equation. Our first goal is 
to prove the same existence and regularity result as obtained by Schoen and Yau [IT], for arbitrary 
positive warping factor (j). This is fairly straight forward, in that the methods of Schoen and Yau 
still apply with little modification. The primary difference is that (especially for applications) it is 
important to track precisely how the estimates depend on <p an d its derivatives. The estimates of 
this section will also be helpful for the case in which <f> is allowed to vanish, but for now we assume 
that (j) is strictly positive. 

Let us set the notation. Consider a hypersurface E C (M xl,j| <p 2 dt 2 ) given by the graph of a 
function t = f(x). This surface has induced metric ~g = g + <j) 2 df 2 , and inverse matrix 

y y 1 + 2 |V/| 2 ' 

where f % = g^ fj and fi = d x if, with x l , i = 1,2, 3, local coordinates on M. Throughout the paper, a 
bar will be placed over geometric quantities associated with E. For instance V will denote the induced 
connection on E, whereas V will denote the connection on the ambient manifold (M xl,j + (f) 2 dt 2 ). 
The unit normal to E is given by 

N= vf-r 2 d t 

v>- 2 + |v/l 2 ' 

and Xi = d x i + fidt, i = 1,2,3 form a basis for the tangent space. Moreover a calculation ([3], [5]) 
shows that the second fundamental form of E is given by 

(2.1) ^sl±M±M±£rhMi 

' " V 1 + * 2 iv/| 2 

In this section it will be assumed that S satisfies an equation of the form 

(2.2) H - TrK = F, 

for some function F to be specified, where H = Ij 13 Aij is mean curvature and TrK denotes the trace 



of K over S. In what follows, as in (|2.2p . we will drop the subscript S when denoting the mean 
curvature and trace operations with respect to E. 

We first estimate the mean curvature and its derivatives. 

Lemma 2.1. Suppose that the surface E satisfies (|2.2p . Then 

(2.3) \H\<\F\ + c{l + \Vlog<f>\), 
and 

(2.4) \VH\ < c(l + |VF| + |Vlog0| 2 + (T^V 2 ^! + 1^1 + 1^1 |V log <f>\), 
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where c is a universal constant. 
Proof. First, by (j2.2|) . we have 

(25) rr..^(, , ftfflWi } F _^ h , PtfflWf? ]F 

( } " 5 [ tJ + yrr^w J " 9 3 a + ^ 2 iv/i 2 ) 3 / 2 + R 

This implies (|2.3[) easily. 

Next, a simple differentiation yields 

1 + ^| V /|2 



=^%;a + d a F + ~Zi h 3 



2# a (^/ Tt )jV/l 2 + 2 (V Q ^)/ ra |V/| 2 + 2 0"V nQ /|V/| 2 
(l + </> 2 |V/| 2 ) 3 / 2 

3 2 (^/ m )|V/| 2 |V7 „ 2 ^ 2 , nv7 
~ 2 (l + 2 |V/| 2 ) 5 /2 Wa|V/l + 2< ^ 7 Vna/) ' 

However, by (|2.ip 

Vy/ = rVl + ^lV/l 2 ^ - - ^-/i - 2 i0'" /,„)/,/• ). 
Thus, after substitution, some cancelations occur and we obtain 

a 9 ^ a+dat + \i+^\vw ^i+www a + 2 iv/i 2 ) 3 / 2 ; 

^ 2 |V/| 2 / W Van^ 2<j>(<t> m fm) 2 fa 0|V/| 2 |V^| 2 /a <I>\V f?{<T fm)<l>a 

(1 + 2 |V/| 2 ) 3 / 2 (l + c/> 2 |V/| 2 ) 3 / 2 (1 + 2 |V/| 2 ) 3 / 2 (l + </» 2 |V/| 2 ) 3 / 2 

; 2^{<rfm)f n , <^|v/i 2 ^ 30 3 |v/i 2 («/» m / m )r 



1 + 2 |V/| 2 1 + 2 |V/| 2 (l + c/> 2 |V/| 2 ) 2 
This implies □ 



We point out that (|2.3j) yields a uniform pointwise estimate of the mean curvature. However, 
(|2.4p illustrates that first derivatives of the mean curvature are estimated in terms of the second 
fundamental form, which is not yet controlled. In fact, estimating the second fundamental form will 
take up a major part of this section. 

We next establish a C° estimate for the second fundamental form of S. This result will arise from 
a Moser iteration applied to the Simons identity. 

Theorem 2.2. Suppose that the surface £ satisfies (j2.2fl . Then 

(2.6) sup \A\ < c 

s 

where c depends on \ log(j)\, | V log 0| , |V 2 log(/>|, and \N(F)\. 

Proof. The proof is lengthy and is divided into several steps. 

Step 1. Derivation of Simons identity and several related inequalities. Recall the Ricci commuta- 
tion formula and Codazzi equations: 

(2.7) V a V n ^4jj V n V a ^4jj — A. m jR ^ an Ai m R j an , 
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and 

(2.8) V a y4ij = VjAy + RNjia- 

Use these to obtain 

AAij = g X7 n X7 a Aij 

= r a (V„Vi Ay + V n R Njia ) 

= 9 (VjVn-Aaj A rn jR ara j A am R j n i + V ' n R]\! jia) 

= (/""^(Vj'^raa + RNajn) ~ d^i^mjR ani A am R j ni — ^ uRn jia) 

= ViVjH + g na (ViRNajn + ^nRNjia — A m jR ani — A am R j ni ). 
Now use the Gauss equations 

(2-9) Rijnd = Rijnd + Ai n Ajd — Ai d Aj n 

to find 

AA i:j =ViVjH + g na (ViR Najn + V n R Njia ) 

9 9 {A m j Rdani + A am Rdjni) 

— 5 <7 Amj^AdnAai — A^A an ) — g na g m ^ A am {A dn Aji — Ai d Aj n ). 

Therefore 

AAij = V l V J H - \A\ 2 A i:j + HA im AJ + G -, 

where 

G-ij = <? (Vj-R/Vajn + V 'nRNjia ~ 9 A m jR da ni ~ 9 A am R d j n i) . 

We then obtain 

\&.\A\ 2 =t l g ia (A da KA ij +g mn V m A da V n A ij ) 

( 2 - 10 ) =g td g ja (g mn V m A da V n A tJ + A^) - | A\ A 

+ g U 9 ]a (HA da A im A] n + A^V^jH). 

This is the Simons identity for S. 
Note that we also have 

iA|A 2 = |^|A|A + |V|A|| 2 , 

so that 

\A\A\A\ = \VA\ 2 - \V\A\\ 2 - \A\ 4 + ^(V.VjH + G l3 + ffi !m A™). 

Set 

(2.11) T :=\VA\ 2 -\V\A\\ 2 . 

For the remainder of this paragraph we will work in an orthonormal basis. Then 

T = £(V n ^-) 2 - \A\- 2 AijVnAj) 2 , 

i,j,n n i,j 

SO that 



|^4| T — — ^ ^ (AijX/ ' n A dm Adm'VnAij) . 



2 

i,j,n,d,m 
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By setting n = i and m = j and using the Schwarz inequality, we have 

|A| 2 T = — '^^(Aij'V iAdj — Adj'ViAij) 2 + — (Aij'V n Ad m — A^ m V n Aij) 2 



2 / j v v 6 My u,j ■ * "J / 2 

i,j,d i,j,n,d,m 



d V i,j IJ / i,j,n,d,m 

Moreover, according to the Codazzi equations (|2.8|) . 

^ ] Aij\7iAdj = ^ ] AijX7dAij + ^ ^ AijRNjdi 

and 

E AdjViAij = AdjVjH + AdjRmji- 

It follows that 

~ 18 E I E A iJ R Njdi ~ E A ^0 H ~ E A dj R Niji J 



2 

i.j.n,d,m 
(n,m)/(i,j) 



after using (a + b) 2 > \a 2 - b 2 . The definition of T in ([2~TT|) now yields 
T>^]T(V„^) 2 



37 

2 



36 1 

w'lX^^E E AijRNjdi - ^2 A dj^3 H - E A dj R Niji 



37 18 



Hence 



+ ^ " 2 { A ij^J nAdm ~ A dmSn A ijf '■ 

i,j,n,d,m 

\A\A\A\ >^ ]T(V n ^) 2 - |^| 4 + ( VaH + G 4j + ^iL4 im A im J 

i,j,n i,j \ m / 



- ^i^i" 2 E E ^i^i* - E a^jj - E a ^ Rn 

d \ i,j j i,j 

2 

+ gyl^l ^ ] (AijV n Ad m ~ AdrnV n Aij) . 

i,j,n,d,m 



2 

\ 
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We can also write 

(2.12) \A\A\A\ > l-Y.^nAj) 2 - \A\ 4 - \H\\A\ 3 + J2^j^jH + 
where 

G\ := Ai J (, ''j ^= • 1 1 X! I Yl A ij R Njdi ~ A dj^jH - A dj R Nijl 

i,j d \ i,j j i,j 

2 

"I - 7^y|^4| ^ ] {Aij\7 n Adm ~ A.(i m \7 n Aij) 



Step 2. An L 4 - estimate for the second fundamental form. Taking two traces of the Gauss equations 
(|2.9p produces 

R = g m g^R l]nrn + H 2 -\A\ 2 . 
Moreover, equation (|2.2p implies that 

R + \A\ 2 = g m g jm R ijnm + (TrK + F) 2 . 

Recall the formula for the scalar curvature of surface £ (see [1], [5]) 

R =2U - J(w)) + \A- K\ 2 + 2\q\ 2 - -divUq) 

+ H 2 - {TrK) 2 + 2K(N, N)(H - TrK) + 2N(H - TrK). 
By combining the previous two formulas we then have 

\A\ 2 + \A- K\ 2 + 2(ji - J(w)) + 2\q\ 2 

(2.13) =-dlv(<t)q) - (2TrK + F)F - 2K(N, N)F 

- 2N(F) + g in g jm R ijnm + (TrK + F) 2 . 

Note that 

\A - K\ 2 = \A\ 2 + \K\ 2 - 2(A,K), 

and set 

G 2 := 2{A, K) - \K\ 2 - 2K(N, N)F - 2N(F) + (TrK) 2 - 2(/i - J(w)) + g in g jm R i3 nm. 
Then, (|2.13p becomes 

2 

(2.14) 2\A\ 2 = -divUq) - 2\q\ 2 + G 2 . 

4> 

Let ip £ C£°(£). Multiply (|2.14p by 4"4> 2 and integrating by parts to find 

/ 2(j)^ 2 \A\ 2 = j -40V(W,9) - 20V 2 k| 2 + <P^ 2 G : 
Jt, Jt, 

' |2 , i„/.2, 



2 



< j 20iwr + ^ g 2 . 

Now replace ip with |j4| p ?/> to get 

(2.15) j UM\ 2+2p < I <t^W(\A\ p n 2 + \U 2 G 2 \A\ 2 K 
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Expand and integrate the first term by parts: 

0|V(L4PV)| 2 = / ^\A\ 2p \V^\ 2 + 2p ( j)^\A\ 2p - 1 ty\A\,V^) + U 2 W\A\ P \ 2 

(j)\A\ 2p \V^\ 2 - -^ 2 A|,4| 2p + ^ 2 |V|^| P | 2 - V> 2 (W>, V\A\ 2p ) 
s 2 



Since 



we then have 



(j)\A\ 2p \Vij\ 2 - W 2 \A\p~A\A\p - 2p^ 2 \A\ 2 P~ 1 ^<p,V\A\). 
|^| p A|,4| p = p\A\ 2p ~ l ~K\A\ + p(p - l)\A\ 2p ~ 2 \V\A\\ 2 , 

<fnly i (\A\ 2 *'- 1 (p'E\A\ + |^| 3 ) + p(p - 1)|A| 2 ^ 2 |V|^|| 2 ) 

< / 0|^| 2p |V^| 2 -2p^ 2 |^| 2p - 1 (V0,V|^|) + i^ 2 |^| 2p G 2 
Jt ^ 

= <P\A\ 2p \V^\ 2 + 2Vl^| 2p (W, V0> + #V| 2p Qg 2 + ^ _1 A^ 

< jf 20i ai 2 ^i wi 2 + <w 2 i4 2p Qc 2 + <r 2 i wi 2 + 0- 1 aV) . 

By using (I2.12p for |A|(A|j4| + \ Af") and setting p = 1 in the above expression, it follows that 

£ ^0V 2 |VA| 2 < jf 2^|^| 2 |V^| 2 + <?tyVl 2 Qc 2 + r 2 |W| 2 + r 1 ^) 

+ / ^ 2 (\H\\A\ 3 -G 1 -A i ^V ij H). 
Jt 

Integrating by parts and applying the Schwarz inequality then yields 

c-V^ 2 |V^| 2 < c4>\A\ 2 \Vi;\ 2 + </>V>V| 2 Qc 2 + c|Vlog</>| 2 ) 
+ / ^ 2 (l^ll^4| 3 -Gi + c|VF| 2 ). 

Jt 

Above and in what follows, c > will always denote an appropriately large constant. Although 
there are derivatives of the Riemann tensor contained within the expression for G%, these may be 
integrated by parts so that 

lj.„/.2|V7 A\2 ^ I A\ /1|2|V7„/.|2 

S 



(2.16) 



Jt 

+ J <p^ 2 (\Riem\ 2 + \VH\ 2 + \V log (j)\ 2 \Riem\ 2 ) 

+ f <pip 2 \A\ 2 (l + \Riem\ + \Vlog<f)\ 2 + \N(F)\ + \V\og(f)\\F\ 
Jt, 

#Vl 3 (l^l + |Vlog^|). 
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Now observe that according to ()2.15p with p = 1, 

HM\ A < I <P\V(\A\ip)\ 2 + ^ 2 G 2 \A\ 2 



< 



jf 20|A| 2 |V^| 2 + [2\V\A\\ 2 + Ig 2 \A\ 2 ^ . 



Also 



Therefore 



|V|A|| 2 = g ij di\A\dj\A\ 

= r(\A\~ 1 A ab V l A ab )(\Ar 1 A cd V J A cd ) 
< c\VA\ 2 . 



or-|.l| J i / 20|^| 2 |VV| 2 + ^V 2 G 2 |^| 2 + c^ 2 |V^P 



Combining this with (|2.16p then yields 



+ y <^ 2 (|iftem| 2 + |VF| 2 + |Vlog^| 2 |ffiem| 2 ) 

+ / (pi) 2 \A\ 2 {I + \Riem\ + \V log cj)\ 2 + \N{F)\ + \V log (f)\\F\) 

J-E 

+ J </>VV| 3 (|#| + |Vlog0|). 
Now replace ^ with ip 2 and apply the Schwarz inequality to obtain 

(2.17) / c -1 # 4 |A| 4 < / ^|V^| 4 + # 4 (l + |iI| 4 + |ViI| 2 + |ffiem| 4 + |Vlog0| 4 + |F| 4 + |iV(-F)| 2 ). 

Let i?^ Q (xo) be a geodesic ball in (M xR,j + cj) 2 dt 2 ) centered at a point xo £ S, with vq less than 
the injectivity radius. Let r be the distance function inMxK and choose the cut-off function ip to 
be a function of r such that if)(r) = 1 for r < ro/2 and ip(r) = for r > tq. Since |Vr| = 1 and hence 
|Vr| < 1, we may further choose if) so that < 1 and < Sr^ 1 . In order to bound the volume 
of S n Bf Q (x ), note that by 

d»Mxi(iV) = # < |F| + c(l + |Vlog0|). 
Upon an integration by parts over the region B^ q {xq) n {(x,i) | i < /(x)}, we obtain 

(2.18) Vol(EnBf (x ))<cr 3 

where c depends on \F\ and |Vlog^|. Applying this result together with ()2. 17[) then yields 



(2.19) / 



\A\ A < c, 



snB 4 rQ (x ) 

where c depends on \N(F)\, \H\, \VH\, \Vlog(f>\, \Riem\, ™*t ' and r^ 1 . This is the desired 
L 4 -estimate for the second fundamental form. 

Step 3. A pointwise estimate of the second fundamental form. Set u = |A| 2 + 1. Then (|2.10j) yields 

An > -c(l + \A\ 2 + \F\ 2 + | V log 0| 2 )u + 2| VA| 2 + 2A ij (VijH + G fj ). 
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Multiply both sides by a nonnegative function £ G C£°(£ n Bi sl (xo)) and integrate by parts to find 

2 

> J (V£, Vu) - c(l + |^| 2 + |F| 2 + |Vlog^| 2 )£u - 2A ij V£V j H 

+ / 2£(| VA| 2 + ^Gij — ViA ij VjH). 

Although Gjj contains first derivatives of the Riemann tensor, these derivatives may be integrated 
by parts. The first derivatives of A which result from this process may be absorbed into |V^4| 2 to 
yield a more simple expression 



0> J V£-Vu + D\Vii)u + D£u, 



where 

D l = -2u- 1 A^V i i? - 2u- 1 A i ig nm RNn jm - 2vr l A nm gV R Nnmj , 

and 

D = -c(l + \A\ 2 + \F\ 2 + |Vlog0| 2 + u- 1 \VH\ 2 + u~ l \A\ 2 \Riem\ + u'^Riem^) 
By ([23D, we have 

il2 , / n 2 



(2.20) sup |£>T + / jD < c. 

EnBt, (x ) JSns 4 ro (jbo) 

We now show that the Sobolev inequality holds on En5^ (xo). The volume estimate (|2.18p allows 
an application of the Hoffman and Spruck inequality [10J for a sufficiently small tq. In particular, 
for any function i/j S C£°(£ n (xo)), we have 

1/3 



QV) <c jj\VM 2 + H 2 ^ 2 ). 



Moreover, by (j2.3|) and Holder's inequality 

c- 1 ( jf V 6 ) < jT | W>| 2 + r 2 (l + |F| 2 + | V log </>| 2 ) Qf V3 . 
Thus, if ro is sufficiently small we obtain the desired Sobolev inequality 

1 /3 

(2.21) (7/) <cjj?M 2 . 

A Moser iteration [S] can now be applied with the help of (12301) and (j2T2"T|) . since 2 > ± dimS = |. 
It follows that 

1/2 



.2 



w(xo) < C ( / It 
snB^, (x ) 

Now the desired estimate (|2.6I) follows from the definition of u and (12.19P . □ 

As a corollary of Theorem 12.21 we have a pointwise gradient estimate for the mean curvature, 
which follows easily from (j2.4p and (j2,6p . 
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Corollary 2.3. Suppose that the surface E satisfies (12. 2|) , Then 

(2.22) sup \VH\ < c 

s 

where c depends on |log0|, | V log 0| , |V 2 log 0| , \F\, and \N(F)\. 

Let y , y 2 , y 3 , y 4 be normal coordinates for the warped product metric near xq £ £, that is 

g = g + 4> 2 dt 2 = g ab dy a dy b , g a b(0) = $ab, d yC g ab (0) = 0. 

Suppose that d y i S T Xo T,, i = 1, 2, 3 so that near xo, £ is given by a graph y 4 = u>(y), y = (y , y 2 , y 3 ). 
Equation (|2.2j) can now be written as 

^ b w a w b \ fv a v b w 



a,b=l 



where W(Y) = w(y) — y 4 , Y = (y 1 ,y 2 ,y 3 ,y 4 ). To see this first observe that Yi = d y i + Wjd y 4, 
i = 1, 2, 3, form a basis for the tangent space to £. Let iV = N a dy a be the unit normal written as a 
1-form. Then 

3 

= g(Y h N) = Y?N a = % N a + w ti N A = N t + W)i N A . 

0=1 

It follows that 

N = El=i^ady a -dy 4 N_ 
\VW\ ' \N\' 

The second fundamental form, in these coordinates, is then given by 

ViNj _ VjVjW 



\N\ \VW\ 

Moreover, the induced metric on £ is g ab — N a N b , with inverse 

K U U W a W b 

-3* - N a N b = g ab - ,„ . 

|VV^| 2 

This yields the desired expression for the generalized Jang equation. 

Parametric estimates may now be obtained as another application of Theorem 12.21 



Theorem 2.4. Suppose that the surface £ satisfies (|2.2p . Then there exist constants p > and c 
independent o/E, ou£ dependent on |log0|c<2,c< and IF^i^, smc/i i/za£ £nl? 4 (xo) C {Y | y 4 = w(y)} 
and | io | £3, a (3 ) < c - 

Proof. Observe that (I2.23D is a strictly elliptic equation for w having the following structure 

3 

E BV {y,w,dw)w y iyj = C(y,w,dw), 
where B t3 (0,0,0) = S 13 . Since the second fundamental form is uniformly bounded and 

|2 ^ f^c w a w c \ f^ d w b w d \ fv a v b w\ fv c v d w 



i^i 2 = E 



9 ~ TvvtTTiT 9 



,vw\ 2 J V lvvri 2 M ivwi M ivwi 

a,b,c,d=l 
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we obtain 

3 

near y = 0. Now by some simple calculus, 

sup (My) | + \dw(y)\ + |«9 2 w(y)|) < c. 

|y|<p 

Therefore the C 3 ' Q estimates follow from Schauder's theory and the C 2 estimates above. 

Lastly, we note that £ may be expressed as a graph for sufficiently small (but uniform) p > 0. This 
follows directly from the pointwise bound on the second fundamental form established in Theorem 
IP □ 

3. The Harnack Inequality 

In this section it will be shown that the quantity {dt , N) satisfies a homogeneous elliptic equation 
with bounded coefficients, when \A\ is pointwise bounded. From this we immediately obtain the 
Harnack inequality as in [T7j, which is used extensively. The proof here consists of a long calculation. 
Although the resulting equation has the same structure as that of [T7], we cannot simply cite this 
reference, as the calculations carried out here must be done in the warped product setting. 

We begin by recalling that the vector fields X{ = <9« + fidt form a basis for the tangent space to 
the surface £ = {t = f(x)}, and that the vector field 

N _ f™ 9 ™ ~ r 2 d t 

vV" 2 + iv/l 2 

is its unit normal. Then 

(dt,N) ' 



V^ 2 + |v/| 2 

Our goal will be to calculate the Laplacian 

A(d t ,N)=r j Vx p Vx ] (dt,N) 

in terms of (dt,N). As usual, in this notation, an over line bar indicates that the particular geometric 
quantity is with respect to the induced metric on S. 

The following lemma may be obtained from the Jacobi equation, as indicated in [2]. However here, 
we carry out a different (and longer) proof, as some of the calculations can be used later. 

Lemma 3.1. For any S C (M xl,g + cj) 2 dt 2 ), not necessarily satisfying (|2.2[) . 

(3.1) A{d t , N) + (\A\ 2 + N(H) + R NN ){d t , N) = 0. 

Proof. First observe that 

A tJ = (X t ,V Xj N), 

so that 

V Xj N = A jp gP*X a . 

It follows that 

V Xj (d t ,N) = (V Xj d t ,N) + (dt, X m )g mn A jn , 

and 

Vx p Vx, (8t , N) = V Xp (V X d t ,N) + g mn (%A jn )(d t ,X m )+ g mn A jn V Xp (d t ,X m ). 
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The Codazzi equations (|2.8p yield 



Moreover 

and 
Thus 



VpA,- n = V n Aj P + RNXjXnXp- 
Vx p {dt,X m ) = X p (d t ,X m ) — F pm (d t ,X a ) 

= (v Xp d t ,x m ) + (d t ,v Xp x m ) -T a pm (d t ,x a ), 



^ x p X m — T pm X a — A pm N. 



A(d t ,N) =g pj V Xp Vx 3 (d u N) 

(3.2) =9 p ^ Xp (V Xj d t ,N) - \A\ 2 (d t ,N) 

+ g mn (V n H + R NXn )(d t , X m ) + AP m (V Xp d t ,X m ). 

The right-hand side will be calculated term by term. 
We first claim that 

(3.3) g mn (V n H + R NXn )(d t ,X m ) = -(d t , N)^ 1 A<j> + N(H) + R NN ), 
To prove this, write 

d t = (d t ,N)N + g mn (d t ,X m )X n 

and observe that 

g mn (V n H + R NXn )(d t ,X m ) = Ric(N,d t ) + 3 t H - (d t ,N)(R NN + N(H)). 
The desired result follows since dtH = and 

Ric(N,d t ) = -(dt,N)(j>- x &<t>. 
To see this last assertion, use the fact that 

(3.4) r| 4 = ri = v\ A = o, rf 4 = (log fa, r| 4 = -<^ 1 , 

to calculate 

R^ijk = 0, RAiij = — ffiij^- 
Here 4 indicates the t-coordinate. This implies that 

Ra = 0, i? 44 = -4>A<f>, 

and hence 

Ric(N,d t ) = -(d u N)Ric(f m d m - ( l ) - 2 dud t ) = -{duN^Acfr. 

This finishes the proof of ()3.3p . 

Next, with the help of (|3.4p . a straightforward calculation yields 

Vx p dt = (log(p) p d t - (t)f p 4> m d m . 

Therefore 

(^X p t ,X m ) = (j)(f m (j) p - f p <t>m). 

Since this is an antisymmetric tensor, and A pm is symmetric, we obtain 

(3.5) AP m (V Xp d t ,X m ) = 0. 

Alternatively, since dt is a Killing field, it follows immediately that (Vx 9j, X m ) is antisymmetric. 
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Lastly, we claim that 
(3.6) 9 pj V Xp (V Xj d t , N) = (d t , N^Acf). 

First observe that since dt is a Killing field, 

(V Xj d t ,N) = -(X7 N d t ,X j ). 

Furthermore 

V N d t = V dt N-L dt N = V dt N, 
where L denotes Lie differentiation. It follows that 

{V Xj d t ,N) = -(V Bt N,X j ). 

Now calculate 

v 9t N = (r 2 + \vf\ 2 r l/2 {fT^da - r 2 n 4 d a ) 

= -(d t ,N}(f i (log<t>) i d t + <t>- 1 <f>%) 
= -(d t ,N)(log<t>) i X i , 

to find 

tf*{V dt N,X j ) = -{d t ,N)Qog<f>y>. 

We then have 

Vx p (r (Vd^Xj)) = X p (gP j (V dt N,X 3 )) +T P pn g^(V dt N,X j ) 



pnz 

(d p (d t ,N))(io g ^r - (d t ,N)d p (\og<t>y - (a t ,iv)r^(iog0)" 

(dt,N)A(log<t>) + (dt,N)(T* -rL)(log0r - (d p (d t ,N))(log<j>y. 



pn pn 

According to a calculation in [4 J (page 760), 



f n A ■ 

-.n .-f n .-S6 n f f 7 ^ 



Moreover, by (|2.ip . 

-d p (d t ,N) = d p (<p- 2 + \Vf\ 2 )- 1 / 2 

(0-2 + | V /|2)3/2 



i" 3 ^ /M pn / n ((log^ P /n + (log</))„/ p + 0/ m m / p /n) 



+ 



Therefore 



2 + | V/ |2)3/2 0-2 + |V/|2 (0-2 + | V /|2)3/2 



(9 t ,iV)(r^ - r^)(log0) n - (9 p (9 i ,iV))(log0y 
(log0)"#™/m/« (log0)^- 3 P 



+ 



(0-2 + |V/| 2 )V2 (0-2 + |V/| 2 ) 3 / 2 

(bg (jiff n [(\Qg 4>)pfn + (log (ft)n/p + M^mfpfn] 
(0-2 + |V/| 2 )3/2 

|Vlog0| 2 



(0-2 + | V/ |2)l/2' 

and (|3.6p follows. 

The desired identity is now obtained by substituting (|3.3p . (|3.5p . and (|3.6p . into (|3.2p . 
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We are now ready to prove the Harnack inequality. 



Theorem 3.2. Suppose that the surface E satisfies (12. 2p , Then there exist constants p > and c 
independent o/E, but dependent on | log (f>\c 2 ,<* & n d \F\ c i, a , such that 

(3.7) sup (d t ,-N)<c inf (dt, —N) 
and 

(3.8) sup |Vlog(«%,-iV}| < c. 

snBf,(x ) 

Proof. The estimates of Theorem 12.41 guarantee that equation (|3.ip is uniformly elliptic for small p. 
Next, we note that the coefficients of (13. 1|) are bounded. To see this, observe that Rnn poses no 
problem, and \A\ 2 is bounded by Theorem 12.21 Moreover 

N(H) 



so that N{H) is bounded by Corollary 12.31 

The Harnack inequality (|3.7p now follows immediately. Standard elliptic theory [9j also guarantees 
that 



sup|V<a t ,JV)| <csup|(5t,JV)|. 



Bp 
1 



Combining this with (|3.7j) yields 



S u V \V{d u N)\<cmf\(d u N)\<cmf\{d u N)\, 

Bp B p Bp 

from which (13.81) follows. □ 



4. Global C 1 Bounds and Existence for the Regularized Equation 

Suppose that the boundary of M is an outermost future apparent horizon with one component, 
that is 

(4.1) e+(8M) = 

where the null expansion 9 + was defined in (|1.7p . All arguments to follow may be easily extended 
to the general case in which the boundary is an outermost apparent horizon, with several future and 
past horizon components. Take a sufficiently large coordinate sphere dooM in the asymptotically flat 
end such that 

(4.2) e + (dooM) > o. 

Following [lj and [15] we construct an extension (M,g,k), of the initial data (M,g,k), having the 
following properties: 

i) M C M with g\M = g, U\m = k, and d^M = dooM, 

ii) g is smooth across dM and k is C ' across dM, 

iii) the region M — M is foliated by surfaces S a , a G [0, do), with 6 + (S a ) < 0, 5o = dM, and 
S ao =dM, 

iv) for all sufficiently small a, 6 + (S a ) < — t? < and k = 0. 
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For convenience, in what follows, we will denote g by g and k by k. Moreover, we will extend the 
warping factor (ft to be positive on all of M. These constructions allow for the existence of appropriate 
barriers, which in turn leads to a solution for the following Dirichlet problem 

(4.3) H(f £ )-K(f £ ) = ef £ on M, 

(4.4) f e = £- on dM, f £ = on djd. 

2e 

The first step is to obtain global C 1 bounds. Observe that a direct application of the maximum 
principle yields the C° bound 

(4.5) sup | / e | < ce" 1 , 

M 

where the constant c depends on \k\ c o and i?. In order to apply the maximum principle to obtain 
bounds on first derivatives, we need to establish the boundary gradient estimates. Let t(x) = 
dist(x, dM), and denote by S T the level sets of the geodesic flow d T = n emanating from dM, where 
n is the unit outer normal (pointing towards spatial infinity) of the surfaces S T . The barrier functions 
ift±(r) will be functions of r alone. In order to find sub and super solutions notice that 

^ (fPifti'ifti \ f (ftVijift + (ftiiftj + (ftjtfti 



'-V 2 A ( (ft^" + 2(ft'^\ , ,_, ^ (ftv tj ift 



+ {g\s T r- 



and 



i + (ft 2 ip> 2 J y^/TT^ft^ J V i + tfip' 2 ' 

\ 9 i + 2 iwi 2 y y " v i + / 

where V'' = • Moreover 

where V 5t denotes the induced connection on 5 r and A St is its second fundamental form. Therefore 
H(ip) - K(ip) - e?ft 



(4.6) (l + «' 2 ) 3 ' 2 v / l + *V" 



(1 + 02^,2)3/2 0Tjr^2 1 + 02^/2 ^ ^l + ^V 2 , 

Let (r) = a — 6r where a and 6 will be chosen appropriately and positive. In order that ift- agree 
with f £ at dM we set a = ^. Then for r € [0, to], with to sufficiently small, 

#(Y>_) - K(ift_) -eift->$-ea + ebr + 0(6 _1 ) > ^ + O^ 1 ) > 0, 

if b is chosen large enough depending on (ft. Similarly an upper barrier may be constructed in the 
form ift+(r) = a + br. Note that in the case of an upper barrier we use the fact that k = near dM to 
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deal with the Trs T k term. Now choose b so large that a + 6ro > sup^y | f £ | and a — bro < — sup^j | f e | , 
then by a standard comparison argument tp_ < f £ < for r G [0,ro]. It follows that 

(4.7) sup|V/ e | < b < ce _1 . 

Moreover, a standard barrier construction at docM also yields a boundary gradient estimate. 

We now use the Bernstein method to obtain global C 1 estimates. For convenience we temporarily 
drop the subscript e from f £ . Differentiate equation (|4.3p with respect to d p to find 



9 lJ 



,J pf;ij "I" 



+ 0;ip/j + <f>if'Jp + 4>:jpfi + 4>jf;ip 



(4.8) 



g 13 



+ <t>ifj + <f>jfi)(<M>p\Vf\ 2 + ^f;m P f m ) 



(l + 2 |V/| 2 ) 3 /2 



— h- ■ 

r"ij;p 



- 2 



tfp- 



f;ij + <t>ifj + 



v./''./'' + o 2 ./;;/' <mw P |v/| 2 + 2 / ;mp / m ) 



i + <^ 2 |v/l' 



(l + ^|V/P) 2 



Notice that the Ricci commutation formula yields 

g'of'L,, = g ij <t>Ptf;pij + R m l3 pf m ) 

= Vj(9 ij <t>f P f;pi)-9 ij <t>f/f;p j 



f P fmR m 



ijp- 



Thus, if u = |V/| 2 , then (gSJ) implies that 



+ + ^u 1 ^ + 



1/2 rinuj^f/f^ 



> £U, 



2Vl + ^ 2 |V/|2 ^ Vl + «^ 2 |V/| 2 

for some coefficients -B* and I?, with bounded. Moreover at a critical point for it, 



9 ij U P 0pf;ij 



'f;pi)=9 i3 U P <t>ph 



7f; PJ )<--cP\V 2 f\ 2 + CU. 



An alternate approach to deal with the g iJ f-ij term is to solve for it in (|4.3h . leaving a manageable 
number of first derivatives of /. Hence from the maximum principle we obtain 

(4.9) sup|V/ e | <ce-\ 



M 



With the global C 1 bounds the equation f)4. 3[) is uniformly elliptic. Standard theory [9] then gives 
C 1,a estimates up to the boundary. The Schauder estimates may now be applied to obtain global 
C 2 ' a bounds. Thus we may apply the continuity method, as in pQ, to the family of equations 

(4.10) H(f g>e ) - 6 K(f gj£ ) = ef Q , £ , 



(4.11) 



u.e 



2e 



on dM, I 



IKE 



< Q < 1, 

on SnoM, 



to obtain a solution of boundary value problem (|4.3p . (|4.4p . Furthermore, by sending dooM to infinity 
we obtain a solution on all of M with the usual decay (|1.9p . Lastly, note that the global estimates 
of this section also hold when M does not have boundary; of course they are easier to prove in this 
case, as there is no need for an extension M or boundary gradient estimates. Thus, in the case that 
M does not have boundary, we may set F = ef e and apply the results of Sections [2] and [3l as well as 
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the global estimates of this section, and let e — > to obtain the desired solution of the generalized 
Jang equation as described in Theorem ll.il (see |17| for details). 

5. Blow-Up Solutions for the Generalized Jang Equation 

Consider the solutions f e given in the previous section, defined on M and with fall-off (jl.9p at 
spatial infinity. Our goal in the current section is to produce solutions which blow-up at the outermost 
apparent horizon boundary of M, by letting e — > 0. We will also construct appropriate super solutions 
in order to obtain an estimate for the rate of blow-up. As in the previous section we will assume here 
that the boundary of M is an outermost future apparent horizon with one component. All arguments 
to follow may be easily extended to the general case in which the boundary is an outermost apparent 
horizon, with several future and past horizon components. 

Observe that the gradient estimate (|4.9p and the boundary condition (|4.4p imply that there exists 
k independent of e, such that 

(5.1) f £ (x) > — for dist(x,9M) < k. 

As in [17J a subsequence of the e-Jang surfaces converges to a properly embedded submanifold of 
MxM. This convergence determines three types of domains: 

f £i {x) — > +oo locally uniformly as 14 oo}, 

f Ei (x) —> — oo locally uniformly as i-> oo}, 

limsup|/ £i (x)| < oo}. 

i— >oo 

By (|5,ip . M_|_ ^ and M + contains a neighborhood of dM. Since dM + — dM consists of apparent 
horizons and the region M — M is foliated by surfaces with 6+ < 0, we must have that M — M C M + . 
Thus dM + is an apparent horizon in M. As dM is the outermost apparent horizon in M, we conclude 
that the closure of M+ is M — M. A standard barrier argument [TTj at spatial infinity shows that the 
e-Jang surfaces are uniformly bounded there, so that Mq contains a neighborhood of spatial infinity. 
It follows that Mq = M, and the limiting Jang surface E blows- up in the form of a cylinder over 
dM. This line of argument was used by Metzger [15] for the classical Jang equation. 

This result holds for all warping factors which are strictly positive up to and including the 
boundary. We would now like to examine what happens when 4> is allowed to vanish on the boundary. 
From now on, let 

t(x) = dist(x, dM), 

the distance to the boundary, and denote by S T the level sets of the geodesic flow d T = n emanating 
from dM, where n is the unit outer normal (pointing towards spatial infinity) of the surfaces S T . We 
assume that 

(f) = r b (f) near dM, 

where <f> is some positive function, and then set 

cf) S = ( T + 5) h 4> near dM, 

where 5 > is a constant. The function cj)g is then naturally extended to the rest of M. When <p is 
replaced by <fis, we write the corresponding generalized Jang equation as 



M+ = {x G M 
(5.2) M_={xeM 
M = {x G M 



(5.3) 



H s (f) - K 5 (f) = 0. 
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For each 5 > a blow-up solution fg of (I5.3P exists. We aim to show that for some subsequence 
5i — > 0, the surfaces X^. = {t = /^(x)} converge smoothly away from the boundary to a blow-up 
solution £ = {t = f(x)} whose asymptotics at the boundary depend on the rates at which 4> and 9 + 
vanish. 

The next result generalizes Theorem 11.21 an d hence a proof of Theorem 11.21 may be found within 
the following arguments. The only difference is that the lower bound is missing. This however will 
be established in the next section, as a by product of the proof of Theorem 11.31 

Proposition 5.1. Suppose that (M,g,k) is a smooth, asymptotically flat initial data set, with out- 
ermost apparent horizon boundary dM consisting of a single future apparent horizon component. 
Suppose further that <ft is a smooth function, strictly positive away from dM and satisfying (jl.8p 
and (|1.12|) . and that satisfies (|1.10|) . Then for b > — there exists a smooth solution f of the 
generalized Jang equation (jl.6p . satisfying (|1.9p . and in a neighborhood of dM 

/<ar- b - i i i + /3 if b>- l —±, 

^ I - 1 

/<-alogr + /3 if b = — , 

for some positive constants a and j3. 

Proof. Let f$ }E denote solutions of the e-regularized generalized Jang equation with warping factor 
<j>S, as constructed in the previous section, 

H S (f S ,e) ~ K S (fs,e) = sf S ,e- 

We now proceed to construct an appropriate upper barrier function ip for fg jE which is independent 
of 5 and e. For this we will assume that (jl.lOp holds for r £ [0, To], with tq sufficiently small. The 
upper barrier ip will be a function of r alone, namely ip = ip(r). Recall that 

„ m ~ m , foV' + WsV , M'0+(S T ) k nn 
H & (vb) - K s {ip) -eif>= 2 . + 



(1 + ^/2)3/2 r 1 + M2 



1 + , = I Tr ST k-ei>. 

In the following, tp will be a decreasing function, that is, ip' < 0. Hence, we have 

M" + 2<f>W k n 



HsW) - K 5 {^) -e*l> = -0+(S T ) + 



1 17/1 



(1 + 0^/2)3/2 1 + ^,2 

(5 ' 5) i h St 

^/l + ^2 ( ^ 1 + 2^2_^/ ) 

Consider 

rb{T)=UT~ a + p, 

for some constants a, a, (3 > 0. In order for \<fnp'\ — > oo as r — > 0, we require that a > b — 1. A 
straightforward calculation shows that 

(i + <^ /2 ) 3 / 2 + 1 + 0^/2 - a 2 a 2 y a + + 1 > T 
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for t sufficiently small. By the assumption 6 + > c r , we obtain 



H s {$) - K s (ij;) -ei>< -c^t 1 + 



C 



(a + b+l)r 



2a-2b+l 



Now set a = b-\ — and notice that this automatically satisfies the previous condition a > b — 1 . If 
a > and a is chosen sufficiently large, then ip is a super solution, 



Kg(ift) -eift<0. 



If a = 0, then setting 



V>(r) 



a log r + (3 



yields a similar result. 

In order to obtain the estimate, recall that \fg e \ will be uniformly bounded away from the boundary, 
say at the surface corresponding to r = tq > 0. This follows from the Harnack inequality and 
parametric estimates, as in [17]. We may then choose j3 sufficiently large independent of 5 and e 
such that ip{jo) > fs, e \s T ■ A standard comparison argument then shows that 



As we have shown above, letting e — > yields a blow-up solution E5 to the generalized Jang equation 
for each positive 5. The solution fg must of course also satisfy the asymptotics (I5.6p . This establishes 
Theorem 11.21 by choosing 6 = 0. 

Let us now extract a subsequence 5i — > such that the surfaces Ylg. converge smoothly away from 
the boundary to a solution E = {t = f(x)} of the generalized Jang equation. This may be proved in 
the usual way |17] , making use of the parametric estimates of Section [2] and the Harnack inequality 
of Section [3j Moreover the asymptotics (|5.6p still hold for /. □ 

Since the warping factor (ft vanishes at the boundary, we cannot say, without further analysis, 
precisely how the solution behaves at the boundary. For instance, although the solutions E^ blow- 
up in the form of a cylinder over dM, it may be the case that as 5{ — > the blow-up solutions become 
arbitrarily close to dM x M and eventually (in the limit) coincide or 'stick' to the cylinder at certain 
points. In order to prevent this, certain inequalities should hold between the vanishing rates of (ft 
and 8 + ; this issue is the primary focus of the next section. 

It turns out that the blow-up rates of (|5.4p are not the only asymptotics for the generalized Jang 
equations. The next result provides other possibilities for upper barriers. 

Proposition 5.2. Suppose that (M,g,k) is a smooth, asymptotically flat initial data set, with out- 
ermost apparent horizon boundary dM consisting of a single future apparent horizon component. 
Suppose further that (ft is a smooth function, strictly positive away from dM and satisfying (|1.8j) 
and (|1.12|) . and that 6 + satisfies (jl.lUp , Then for b > 1/2, there exists a smooth solution f of the 
generalized Jang equation (jl.6p . satisfying (|1.9p . and in a neighborhood of dM 



(5.6) 



fs,e\s T < V'(t) for T€(0,T O ]. 



(5.7) 




for some positive constants a and j3. 
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Proof. Recall the calculation in (I5.5p . Our strategy will be to choose the upper barrier function ^s, 
so that the second derivative of iftg will dominate all other terms in (|5,5|) . except — 0+(S T ). To this 
end, consider 

+ = M'S 

T + 

for some A to be determined. It is an easy exercise to show that its solution is given by 

f 1 e Xs 

(5.8) ^(r) = /i 1 / + 

Jt+8 s 

for some positive \i\ and fi2 to be determined. It is obvious that ip$ is monotonically decreasing. 
Therefore 

+ 2(log foYti = if>'i + -^-z + 2(log 4>)'i>' s < tfl + -^tfs - c$ 

T + T + 

= (A - cW & < ^ < 



for large A. Hence, (15. 5j) becomes 



- - e^s < -e + + — , "S 2 , 3/9 + c^V«l 



2(1 + <^ 2 ) 3 / 2 



<-c-v- u-cj ^v«r 2 <o. 

if again A is large enough. It follows that ^ is a super solution. 
Recall that fs, s \dM < ce _1 by (|4.5p . and note that 

^(0)>a^~ 26 + /3 if 6>^, 
^(0) >-alog<5 + /3 if 6=1 

This shows that V^ldM > /<5,e|dAf for sufficiently large a, if 5 = e^ 26-1 ) 1 when 5 > 1/2, and 5 = e~ £ 
when b = 1/2. The maximum principle then implies (if /3 is chosen large enough) that fg jE < tpg for 
r G [0, ro]. Let us now extract a subsequence e, — >• such that the surfaces Ej converge smoothly 
away from the boundary to a solution £ = {i = /(x)} of the generalized Jang equation. This may 

be proved in the usual way as at the end of Proposition 15.11 even though here 5i = ef h ^ when 
b > 1/2, and Si = e~ Ei when b = 1/2. It follows that the solution / must satisfy the asymptotics 

□ 

Again, because the warping factor 4> vanishes at the boundary, we cannot say, without further 
analysis, precisely how the solution behaves at the boundary. It might coincide or 'stick' to the 
cylinder over the boundary at certain points. Nevertheless, the barrier construction yields an upper 
bound for the asymptotics. 

We may now compare the different asymptotics (|5.4p and (|5.7p . Observe that 2b — 1 < b + 
only when b < and therefore the asymptotics in (|5,7p only improve the ones in (j5.4j) when 

\<b<^. 
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6. Further Analysis of the Blow-Up Solutions 

The purpose of this section is to obtain subsolutions for the generalized Jang equation and to 
prove Theorem ll.31 It will be shown that subsolutions exist with the same asymptotics as described 
in Propositions 15.11 and 15.21 An appropriate comparison argument will then be employed to prove 
Theorem ll.31 As in the previous two sections, we will assume for simplicity that the boundary of M is 
an outermost future apparent horizon with one component. All arguments to follow may be extended 
to the general case in which the boundary is an outermost apparent horizon, with several future and 
past horizon components. Moreover the additional arguments needed for such an extension will be 
recorded at the end of this section in the proof of Theorem 11.31 

Let ip(r) be a function of r alone, which satisfies ip' < and 

(6.1) \<jnp'\ -»• oo as r -> 0. 

Then according to (|5.5p there exist bounded functions c\ and C2, with c\ > | for small r, such that 

(6.2) - KM = -0 + (Sr) - ^ + |) + 

Lemma 6.1. (1) Suppose that - l -r± < b < and that (fTTTOl) and ([02]) hold for r G [0,t ]. Then 
for sufficiently small To there exist sub and supersolutions ip and ip of the generalized Jang equation, 
satisfying (|6.ip . and such that 

a- 1 T- b - L r + /3- 1 <^,^<ar~ b -^+/3 if - — < b < 

(6.3) ( _ ?_l 2 
—a 1 logT + /3 1 < ip, if} < — alogr + /3 if b = — , 

for some positive constants a and /3. Moreover 

H{±) - K(f) > \t\ H@) - K@) < -^-t 1 , 

where the constant A > may be chosen arbitrarily large and c _1 is as in (jl.lOp . 

(2) Suppose that \ < b < ^fci, and that (I1.12p and (jl. 13j) hold for r G [0, To]. Then for sufficiently 
small To there exist sub and supersolutions x andx of the generalized Jang equation, satisfying (|6,ip . 
and such that 

a^T 1 - 2b + r 1 <X,X<^ b + /3 if \<b< 1 -^-, 

-a" 1 logT + /3" 1 < x,X < -alogr + /3 if b = -, 

for some positive constants a and f3. Moreover 

H(x) - K(x) > At 7 , H(x) - K(x) < -A/, 

where the constant A > may be chosen arbitrarily large and 2b — 1 < I < min(i, 2b). 

Proof. (1) First consider the case — ^ < b < ^jp and take the most obvious choice 

iP(T) = ar- a + (3, 

where a is to be determined. A calculation shows that 

ip" 2b _ 2b - a - 1 
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and \<fr<p'\ = actT b ~ a ~ 1 — )• oo as r — > if a > b — 1. By applying (jl.lOp to (|6,2p we then have 

H&) - K{#) = -c 3 r l - (2fe-a-l)c V a _ 2fe+1 + ^2 r 2a- 2b+2) 

a z ar arar 

where C3 > c _1 > is a bounded function. Setting a = b + yields 

Thus if To is sufficiently small, then by defining ift to have the structure of 1/1 with a sufficiently small 

H(f)-K(f) > a~V, 
and by defining ijj to have the structure of ip with a sufficiently large 

H($)-K($)<-^-t 1 . 

Now consider the case when b = — and take 

V'(r) = — q log r + /?. 

Observe that 

?//' 26 _ I 

and | 1 = ckt "3~ — > 00 as r — > 0. By applying (jl.lOp to (|6.2p we then have 

ffty) - K(V) = -c 3 r' + + ^|r' +1 . 

or a z 

Thus if ro is sufficiently small, then by defining ip to have the structure of ip with a sufficiently small 

H{f)-K(f) > q-V, 
and by defining -0 to have the structure of tfj with a sufficiently large 

H$)-K$)<- C —T l . 

(2) Assume that ^ < b < and let x( r ) De a function of r alone satisfying (I6.ip . According to 
and ([TT3D 

H(X) " * 00 = c 4 r' - ^ + I) + 
where C4 is a bounded function. This suggests that we study the ODE 

The choice of =fA will be used when defining the sub and supersolutions, respectively. In order to 
solve (]6.4p observe that this equation is equivalent to 

-4b \ ' 

= ±2\t 1 - 2 \ 



X' 2 



and hence 



.5) X = - ~ + A z t 

\l + l-2b 



1 



/ ±2A l+l+2b 1 A 2^4b\ 2 
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for some constant A. We choose A > 0, since if A = and —A is chosen in (16. 4|) . then (|6.5p yields 
similar asymptotics as in (1). Notice also that the expression inside the square root is positive for 
small t, since 46 < I + 1 + 26 and A > 0. It follows that 

^^ 2h < -X' < 2A- 1 r- 2b 
for sufficiently small r. This shows that (|6.ip holds. Furthermore 

H( X ) - K( X ) = c 4 r l ± ci\r T + c 5 A 2 r 2fe , 

for some bounded function C5. We define Xi X to be the solutions of (|6.4p constructed above and 
corresponding to —A, +A respectively. By choosing A sufficiently large the desired result follows, 
since I < min(7, 26). □ 

The existence of two subsolutions with different asymptotics, when i < 6 < -rr, indicates that 
there will be two different blow-up solutions of the generalized Jang equation, one corresponding to 
each of the distinct asymptotics. These solutions will arise from two different sequences of solutions 
to the regularized equation. More precisely, consider the generalized Jang equation with (j) = T b 4> 
replaced by (f)$ = (t + 5) b cf), as in (|5.3p . By Theorem 11.21 for each 5 > 0, there exists a blow-up 
solution t = fs which, according to the proof of Proposition 15.11 asymptotically approaches the 
cylinder dM x K at a rate given by (|5.4p . A subsequence will then converge to a blow-up solution 
of the generalized Jang equation as 5 — > 0, and this solution will satisfy the asymptotics (|1.14p . 
The second sequence of solutions will arise from the (^-regularized generalized Jang equation, and 
will be constructed to have finite values at dM. However, as 5 — > these boundary values will 
become arbitrarily large, and a subsequence will converge to a blow-up solution of the generalized 
Jang equation which satisfies the asymptotics (|1.15p . 

We first analyze the case of the solutions t = f$ with asymptotics given by (|5,4p . Consider the 
(^-regularized generalized Jang equation applied to a function ip(r) as in (j6.2[) . 

( 6 . 6 ) Hl wo - k sW = -e +( s T ) - I7T £L_ + + 

Let be the subsolution constructed in (1) of Lemma IfTl] and consider 

It is clear from the proof of Lemma 16.11 that ^ is a subsolution of the (5-regularized generalized Jang 
equation, 

(6.7) H 5 (± s ) - Ksi^) > A(r + 5) 1 > 0. 

We will show that ip acts as a lower barrier for the regularized solutions fg. 

Proposition 6.2. If fs\s TQ > ^ 5 ( r o) then fs > ± s for all r G [0,r ]. 

Proof. Consider the function w$ = fs — ip s - Since V^(0) is finite and fs blows-up at dM, we have 
that ws\dM > 0; moreover, by assumption ws\s TQ > 0. Suppose that w$ attains an interior negative 
minimum, then at that point 

It follows that at the minimum point 



V/a = V^ 5 , V 2 w 5 > 0. 



(Hs(fs) - Ks(fs)) - (Ha(^) " K S ® 6 )) = U - j^MttA , > 0- 

V 1 + Ji + <pj\vfs\ 2 
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On the other hand, by (16. 7p 

(Hs(fs) - Ks(fs)) - (Hsil,) ~ K ^±&)) ^ " A ( T + < °- 
This contradiction shows that w$ > for all r E [0, To]. □ 

We are now ready to prove the primary existence result for the first asymptotics of Theorem 11.31 

Theorem 6.3. Suppose that (M,g,k) is a smooth, asymptotically flat initial data set, with outermost 
apparent horizon boundary dM consisting of a single future apparent horizon component. Suppose 
further that <f> is a smooth function, strictly positive away from dM and satisfying (j 1 . 8 1) and f j 1 . 1 2 1> . 
and that 9 + satisfies (jl.lOp for r £ [0,ro]. // — b < 2 ^ en ^ iere exists a smooth solution f 
of the generalized Jang equation (jl.6p . satisfying (jl.9p . and such that the following estimates hold in 
a neighborhood of dM 

a- 1 r- 6 -^ i +r 1 </<ar- b -¥ + /3 if <&<— , 

(6.8) 2 2 
-a^logr + r 1 < f < -alogr + /3 z/ 6= — , 

for some positive constants a and f3. 

Proof. By Theorem 11.21 for each 5 > 0, there exists a blow-up solution t = f$ of equation (|1.6p with 
</> replaced by 4>$, and this solution asymptotically approaches the cylinder dM x M at a rate given 
by ([5.4p , As in the proof of Proposition 15.11 there is a subsequence <5j — )■ such that these graphs 
converge smoothly away from the boundary to a solution t = f of the generalized Jang equation 
(II. 6p . which satisfies (II. 9p . This implies that there exists a fixed value for ip (to), independent of <5j, 
such that fsi\s T > V'r ( r o)- I n fact by writing 

we can achieve this by choosing fixed a« = a and j3i = f3 for all i. Therefore ip converges to a 
function ip satisfying the estimates (|6.3j) . Moreover by Proposition 16. 2\ fg i > ip s for r £ [0, to], and 
therefore / satisfies (j6.8j) . □ 

Blow-up solutions of the generalized Jang equation satisfying (|1.15p will now be constructed. 
These will arise from a sequence of solutions to the ^-regularized equation with finite boundary 
values. Assume that ^ < b < and that (j!.13j) holds, and let Xi X De t ne functions constructed in 
Lemma IBTTl Set 

X s (t) ■= X(j + 8) and Xs(. T ) '■= x{ T + <5)> 
and note that these translated functions are sub and supersolutions for the ^-regularized generalized 
Jang equation. In particular, from (|6.6p and the proof of Lemma 16. II we have 

(6.9) H 5 (x s ) - K s (x s ) > A(r + 5) T , H s (x 5 ) - K 5 {x & ) < -A(r + 6) T . 

Let Ts = (5 1_2fe if b > \ and T5 = — log 5 if b = |, and observe that T$ ~ ^(O) ~ X<5(0). Consider 
now the following boundary value problem for the (5-regularized generalized Jang equation 

(6.10) Hs(h s ) - K 5 {h s ) = on M, 

(6.11) hg = Ts on dM, h$ — > as |a;| — > 00. 
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Proposition 6.4. Suppose that dM is an outermost future apparent horizon having one component, 
and that (|1.13|) holds for r G [0,To]. If \ < 6 < ^jp, i/ien i/iere exists a smooth solution of (|6.1U|) . 
(|6,lip . which is C° up to the boundary, and satisfies the following estimates in a neighborhood of 
dM 

a~ 1 (r + 5) 1 - 2b + /T 1 <h 5 <a{T + 5) l - 2b + (3 if I < 6 < i±l, 

(6.12) 2 ] 2 
-a" 1 log(r + 5) + /T 1 < hs < -alog(r + <5) + /? if b = -, 

for some positive constants a and ft independent of 5. 

Proof. We first solve a boundary value problem for the ^^-regularized generalized Jang equation: 

(6.13) H s (hs, E )-K s {h S)e )=e(h s>e - X ) on M, 

(6.14) hg >e = Ts on dM, hs )£ — > as \x\ — >■ oo, 

where x is a smooth function on M such that x = T& on sup Af |x| < and x = in a 

neighborhood of spatial infinity. Moreover, we require that \ is decreasing in r near dM. 

We now employ the continuity method to show that a unique solution exists for (|6.13|) and (|6.14p . 
This is similar to arguments used for (|4.1Up and (|4.11|) . To begin, observe that the maximum principle 
yields the bound 

SUp|/l,5 ie | < T S + CE~ X . 

M 

To establish boundary gradient estimates, we construct subsolutions and supersolutions. 

Fix a ro sufficiently small. For subsolutions, take x Se as in Lemma 16.11 Specifically, x satisfies 

H s (x St£ )-Ks(x St£ )>X(r + S) T , 

for some positive A. Since this subsolution is a solution to a second order ODE, there are two free 
parameters that may be chosen appropriately so that 

X 5i£ (0) = T 5 , x^(t ) = -(T S + ce- 1 ), 

and 

Xg >e <X for t £ [0,r ]. 

This shows that wg tE := h$ )£ — x Se — along the surfaces dM and S TQ . Suppose that ws. £ achieves 
an interior negative minimum. Then at that point, a standard comparison argument yields 

[Hs(h s , £ ) - K 5 (h 5 , £ )] - [Hs(x s J ~ K S (XsJ] > 0- 
On the other hand, according to (|6.9p and w§ e < at the point in question, we have 
[Hs{h s , £ ) - Ks(h 5 , £ )} - [Hs(x 5 J ~ Ks(x 5:£ )} < <h &>£ - x ) - A(r + 5) T 

= ew 5 , £ + e(x St£ ~ X) ~ A(r + 5) 1 < 0. 

It follows that hg j£ >Xg £ f° r au T ^ [0> r oL an d this yields a lower bound for d T hs t£ . 

For supersolutions, let X be defined in the same way as x> from the proof of Lemma 16. H except 
that (]6.5p is replaced by 



i 

-2A 
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Then xs,e '■= x( T + <5) satisfies 

Hs(xs,e) ~ K S (X5,e) < -A(r + S) T 
Choose the two parameters defining xs.e appropriately so that 

X5,e(0) = Ts, XS,e( T o) = T S + C£~ X . 

Then Xg e — X f° r an T £ [0,to]. As above, a comparison argument can be employed to show that 
hs, £ < X<5,£ for all r G [0, ro]. This yields an upper bound for d T hg :£ . A solution of (|6.13p . (|6,14p is 
now guaranteed, and uniqueness follows from a maximum principle argument. 

As in the proof of Theorem 16, 3\ we may now let £ — > and extract a subsequence of solutions, still 
denoted hg >£ , which converges smoothly away from the horizon to a solution hg of (|6. 10p . Moveover 

hg — > as |x| — >• oo. 

We need to show that 

hg = T s on dM. 

Note that the functions hg >e are uniformly bounded, independent of e (in fact independent of 5 as 
well), when restricted to the surface S TQ . It follows that hs j£ \s T0 < Ts for 5 sufficiently small. We 
may then choose sub and supersolutions x & ^Xs-, satisfying (16. 9p . the estimate (16.121) with constants 
a and /3 independent of e, and such that 

X 5 (0) = hs, e | dM = Xs (0) , Xj (t ) < h &)£ \s TO <X5( T o), 

and 

X s < X < Xs for T G [°i r o]- 

A comparison argument, as in the first half of this proof, may now be used to establish (|6.12p for 
hs s . This establishes C° and boundary gradient estimates, independent of e. We may now apply 
Theorem 15.2 of [9] to obtain global C 1 bounds independent of e. It follows that the limit hg is C° 
up to the boundary. 

Lastly we may take the subsequential limit hg — > h as 5 — > 0, to find that the functions hg are 
uniformly bounded, independent of 8, when restricted to the surface S T0 . Then as in the previous 
paragraph, sub and supersolutions x s and Xg m ay be used to obtain the estimate (I6.12p with constants 
a and f3 independent of 5. □ 

Theorem 6.5. Suppose that (M,g,k) is a smooth, asymptotically flat initial data set, with outermost 
apparent horizon boundary dM consisting of a single future apparent horizon component. Suppose 
further that 4> is a smooth function, strictly positive away from dM and satisfying (jl.8p and ()1.12p . 
and that 9+ satisfies (|1.13|) . If ^ < b < then there exists a smooth solution h of the generalized 
Jang equation (|1.6|) . satisfying (|1.9p . and such that the following estimates hold in a neighborhood of 
dM 

a -V- 26 + /T 1 < h < aT l ~ 2b + /3 if - < b < L±A 
(6.15) 2 ! 2 

-a _1 logr + /3 _1 < h< -alogr + /3 if b = -, 

for some positive constants a and j3. 
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Proof. By Proposition 16.41 for each 5 > there exists a solution t = hg of boundary value problem 
(|6.1Up , (|6.11|) . As in the proof of Theorem \ti.'S\ there is a subsequence Si — > such that these graphs 
converge smoothly away from the boundary to a solution t = h of the generalized Jang equation 
(jl.6p . which satisfies (jl.9p . Since the constants a and /3 in (|6.12p are independent of 5, the limit h 
satisfies (I6TT5D . □ 

Proof of Theorem \l.S\ This follows from Theorems 16.31 and 16, 51 if the boundary has one future appar- 
ent horizon component. Similar arguments hold if the boundary has a single past apparent horizon 
component. 

Now consider the case of multiple components in which some components belong to the category 
(1) and some belong to the category (2). The asymptotics of (1) arise from taking the limit as e — > 
to obtain a blow-up solution, and then taking the limit as 5 — > 0. The asymptotics of (2) arise from 
taking the limit as e — > to obtain a solution with finite boundary values ()6.1ip . and then taking 
the limit as 5 — > to obtain the blow-up solution. Moreover the relevant constructions may be 
performed in a neighborhood of each component, so that this process of letting e — > first, and then 
taking the limit 5 — > 0, may be carried out to obtain the desired result. □ 



We now give an example to show that the nice asymptotics described in Theorems 16.31 and 16.51 
can fail if b = ^fc^. In particular, we will exhibit a solution of the generalized Jang equation which 
sticks to the cylinder. 

Example 6.6. Let us consider the exterior region of the Schwarzschild spacetime with metric 

r J \ r J 

where da 2 is the round metric on § 2 . Let t = fir) be a radial graph, with induced metric 

^((i-vr-^-vH* 2 ^ 2 - 



As is calculated in [2], the second fundamental form of the graph is given by 

ki 



ffiijf + cf>ifj + (fijfi 



where <p = y 1 — ^ and the covariant derivatives are calculated with respect to the metric g. Thus 

(M = M 3 — B2m{0),g, k) forms an initial data set for which the graph t = f(r) is a solution of the 
generalized Jang equation. We choose a solution t = f(r) such that the function / is smooth up to 
the boundary r = 2m. In particular, /(2m) < oo and we may consider this as an example of a Jang 
graph 'sticking' to the cylinder. For such an /, we note that 

V r I \ r I V r 



Therefore the distance to the boundary is given by 
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This implies that cj) ~ r, or rather 6 = 1. We now compute the null expansion of the coordinate 
spheres with respect to the initial data metric g. A standard formula [4| yields 



Vgl 1 / 2m 
-tfsv = ~ \ l , 

and the trace of the initial data k over the coordinate spheres is given by 

<h^\f 2 (jig 11 ?' ( A 2m \ § 



Tr Sr A; = ^ = = = =^= = ~ 1 



It follows that 



2 

0± =H ST ±Tr Sr k~ a/1-— • 
In terms of the notation used above, we have #± ~ r, and in particular, 1 = 1. □ 
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